Scaling properties of three-dimensional magnetohydrodynamic turbulence 
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The scaling properties of three-dimensional magnetohydrodynamic turbulence are obtained from 
direct numerical simulations of decaying turbulence using 512 modes. The results indicate that the 
turbulence does not follow the Iroshnikov-Kraichnan phenomenology. The spectrum is consistent 
with k~^'^. In the case of hyperresistivity the structure functions exhibit a clear scaling range 
yielding absolute values of the scaling exponents ^p, in particular (^3 ~ 1, consistent with a recent 
analytical result. The scaling exponents agree with a modified She-Leveque model (^p = p/9 + 1 — 
(1/3)''''^, corresponding to Kolmogorov scaling but sheet-like geometry of the dissipative structures. 
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Magnetic turbulence is the natural state of a plasma 
in motion, especially in astrophysical systems. The con- 
venient framework to describe such turbulence is magne- 
tohydrodynamics (MHD). For high magnetic Reynolds 
number Rm = vIq/t], where u is a typical turbulent ve- 
locity, Iq the integral scale and 77 the magnetic diffusivity, 
there is a broad range of scales I between Iq, and the dis- 
sipative scale length l^, Iq ^ I :^ lii, called the inertial 
range which exhibits characteristic self-similarity or scal- 
ing properties. 

The concept of inertial-range scaling was introduced 
by Kolmogorov ||i| for hydrodynamic turbulence, which 
is called the Kolmogorov (K41) phenomenology. As- 
suming homogeneity and isotropy of the turbulence as 
well as locality of the turbulent cascade process, he ob- 
tains e ~ {5viY /ti = {5viY /I yielding the scahng law 
5vi r^ e^'^l^'^ . Here e is the energy dissipation rate and, 
to be specific, Svi — [v(x + 1) — v(x)] • \/l is the longi- 
tudinal velocity increment. A direct consequence is the 
Kolmogorov energy spectrum Ek '^ e^/'^fc"^/^. 

For MHD turbulence the Iroshnikov-Kraichnan (IK) 
phenomenology g, [^ takes into account the Alfven ef- 
fect, the coupling of small-scale velocity and magnetic 
fluctuations by the integral-scale field Bq. Hence the 
natural variables are the Elsasser fields z* = v ± B, 
which describe Alfven waves. In the IK phenomenol- 
ogy the spectral transfer is reduced by the factor ta/ti, 
e ^ {ta/ti){5ziY /tu where ta = 1/va, va = Alfven ve- 
locity in the field Bq, Szi ~ 5vi '^ SBi ^ (e-w/i)^/''^^/*, and 
the IK energy spectrum becomes Ek ^ {evAY^^k'^^'^. 
This spectrum can also be written in the form Ek ~ 
g2/3^-5/3^^^^-ji/6 .^^j^j^ ^jjg integral scale ^o defined by 

^0 = ^a/^' which illustrates the nonlocal character of 
the energy cascade in MHD turbulence. 

It is, however, well known that these qualitative scaling 
relations for Svi or 6zi are not exactly valid in a statisti- 
cal sense because of intermittency, which implies that the 
distribution of turbulent scales is not strictly self-similar. 
A quantitative measure is provided by the scaling expo- 
nents (p of the structure functions, the moments of the 



field increments. For hydrodynamic turbulence She and 
Leveque B proposed a model leading to the expression 
(|L ^p/Qj^2[l- (2/3)P/^], which fits the experimental 
results surprisingly well, reproducing in particular the 
exact result C3 = 1. This model has been modified for 
MHD incorporating the IK effect g, ^, which yields 
Cp^ = p/8 + 1 - (l/2)P/4, in particular (f^ = 1. 

The IK phenomenology has been supported by direct 
numerical simulations of 2D MHD turbulence at moder- 
ate Reynolds numbers |Q. However, recent developments 
in MHD turbulence theory cast some doubt on the gen- 
eral validity of the IK scaling. 2D simulations at con- 
siderably higher Reynolds numbers reveal an anonralous 
scaling behavior |^ , J9J , indicating that the results of Ref. 
H are not asymptotic. There have also been theoretical 
arguments in favor of a Kolmogorov scaling, e.g., pO| , 
pT| . Even more convincingly, exact relations have been 
derived for moments of certain triple products of Szi |l2| , 
which are shown to be proportional to Z, i.e., Cs = li anal- 
ogous to the well-known 4/5-relation in hydrodynamic 
turbulence, thus excluding the IK result ^4 = 1. Scaling 
exponents for MHD turbulence have also been derived 
from observations in the solar wind |13| . Here agreement 
with the IK exponents has been claimed |g], but in this 
comparison the observational results were normalized as- 
suming ^4 = 1. Actually the error bars seem to be too 
large to reach a definite conclusion. 

To clarify the issue of scaling in 3D MHD turbulence 
direct numerical simulations are desirable with higher 
Reynolds numbers than studied previously, for instance 
in ||I||, [[5), (l|], 0, ||I||. In this Letter we present a nu- 
merical study of freely decaying turbulence with spatial 
resolution of 512'^ modes. The scaling properties are ana- 
lyzed by considering the time-averages of the normalized 
spectra and structure functions. We solve the incom- 
pressible MHD equations 



9tB - V X (v X B) = r]^{-l)''-^V'^''B 



(1) 



9tw - V X (v X w) - V X (j X B) = H^{-lY~'^V"'w 



(2) 



w = Vxv, j = VxB, V-v = V-B = 0, 



by applying a pseudo-spectral method with spherical 
mode truncation as conveniently used in 3D turbulence 
simulations (instead of full dealiasing by the 2/3 rule 
used in most 2D simulations). The generalized magnetic 
Prandtl number rj^/ ^u has been set equal to unity. Initial 
conditions are 
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which are characterized by random phases at, /?k and 
satisfy the conditions k-Bk = k-Vk = 0, £■ = E^+E'^' = 
1 and E^ /E^^ — 1. Further restrictions on Bk arise 
by requiring a specific value of the magnetic helicity 
H = J d'^x A-'B. We believe that finite magnetic helicity 
is more typical than i? ~ 0, since MHD turbulence usu- 
ally occurs in rotating systems. The wavenumber fco, the 
location of the maximum of the initial energy spectrum, 
is chosen fco = 4, which allows the inverse cascade of H^ 
to develop freely during the simulation time of 10 time 
units (about 7 eddy turnover times, defining the eddy 
turnover time as the time required to reach the maxi- 
mum dissipation from the smooth initial state). Though 
this choice implies a certain loss of inertial range, the 
sacrifice is unavoidable in the presence of inverse cascade 
dynamics, since assuming fco ~ 1 would lead to magnetic 
condensation in the lowest-fc state, which would also af- 
fect the dynamics of higher-fc modes. Both normal dif- 
fusion i^ = 1 and hyperdiffusion v = 2 have been used, 
i^ = 1 to discuss the spectral properties and i^ = 2 to de- 
termine the scaling of the structure functions. All runs 
presented in this Letter have finite H, H/Han^^ ~ 0.6, 
and negligible alignment. Table | lists the important 
parameters of the simulation runs, where the magnetic 
Taylor Reynolds number is Rm\ — Rm}'^. Since Rm\ 
is not stricly constant during turbulence decay but in- 
creases slowly Ruix ^ t^/^ , we give the values taken at a 
specific time t = 4. 

We first discuss the spectral properties considering the 
angle-averaged energy spectrum Ek ■ Figure H shows the 
scatter plot of the normalized spectrum (the normaliza- 
tion is discussed below), compensated by fc'^'^, taken 
from run 3 over the period i = 4 — 10 of fully developed 
turbulence. The spectrum exhibits a clear scaling range 
of almost one decade with a spectral law, which is defi- 
nitely steeper than the IK spectrum fc^"^/^, close to (in 
fact slightly steeper than) k^^/^ indicated by the dashed 
line. In order to form the time average in a system of 
decaying turbulence, the spectrum must be normalized 
to eliminate the time variation of the macroscopic quan- 
tities. In hydrodynamic turbulence the only such quan- 
tity is e, which leads to the universal form of the Kol- 
mogorov spectrum E{k) = Ek/{erf'Y/^, k = Mk, where 



Ik = ('7'^/e)^^^ is the Kolmogorov length. However, when 
normalized in this way the MHD energy spectrum is 
found to change during turbulence decay and even more 
strongly so when comparing runs of different Rm\. For 
finite magnetic helicity the spectrum may also depend on 
H, which introduces a second macroscopic length scale 
li — H/v\ in addition to Iq, i.e, the spectrum may con- 
tain some function of lo/li. To determine this function 
we propose the following argument. Since the Alfven 
effect is clearly present in the simulations, kinetic and 
magnetic energy spectrum being nearly equal at small 
scales, while on the other hand the scaling Ek ^ k^^^^ 
is observed, we modify the nonlinear transfer in the IK 
ansatz by a factor {l/liY , 
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and determine 5 by requiring the observed scaling Szi 
l^/"^, which gives 5 = —1/3 and hence 
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Also the dissipation scale length is slightly changed. Bal- 
ancing nonlinear transfer and dissipation gives 



^K — ^K{h/h 



,1/8 



(6) 



Using these relations we obtain a the normalized energy 
spectrum E(k) = Ek/[{eTff'^{lolhf/^] = k-^/^F(k), 
k = kl^ . Normalized in this way the spectra at different 
times of run 3 coincide very well as seen in Fig. ^ and so 
do the time-averaged normalized energy spectra of runs 
1-3 shown in Fig. g, which vary only by the extent of the 
inertial range, apart from statistical oscillations. Rela- 
tions (^, (|6|) are not valid for _ff ~ 0, where we expect 
the pure Kolmogorov normalization to be valid. 

A more complete picture of the inertial-range distri- 
bution of turbulent structures is provided by the scaling 
exponents (^p of the structure functions, where the second 
order exponent ^2 is related to the inertial-range spec- 
tral law r^ k^^^^''^\ To be definite we consider the mo- 
ments of the absolute value of the longitudinal increments 
6zi discussing only the runs 4 and 5 with the highest 
Reynolds numbers. The normalized structure functions 
S^{T) = {\Sz^\P)/EP^^, T= l/l^, are averaged over time, 
exhibiting a similar weak scatter as for the spectrum in 

Fig. |i|. 

For normal diffusion 1/ = 1 no scaling range is visi- 
ble. (Note that the structure function ^2 corresponds to 
the one-dimensional spectrum Ek^ , which has a shorter 
inertial range than the angle-averaged spectrum shown 
in Fig. ^) For 1^ — 2, however, there is a scaling 
range 30 < I < 200, as seen in Fig. S, where the time- 
averaged curves S^ are plotted for p = 3,4. The in- 
serts give the logarithmic derivatives, where the central 



quasi-constant parts determine the scaling coefScients, 
the dashed horizontal lines indicating the most probable 
values (^ ~ 0.95, Ct — 1-15. These results are consistent 
with the spectral law derived from Fig. m and are close 
to the analytical prediction (^3 = 1. It is true that the 
analytical theory refers to third-order moments different 
from (|i52:|^) discussed here, but the scaling coefficients 
should not depend thereof (the scaling range, however, 
does). One might object that the use of hyperdiffusion 
affects the inertial-range scaling, if the scaling range is 
not very broad. In fact, the energy spectrum law tends 
to be polluted by the bottleneck effect, which is particu- 
larly pronounced for hyperdiffusion (see e.g., |J9|). Thus 
the energy spectrum in run 4 (not shown) is effectively 
flatter than fc"^/'^ expected from the value of (2- How- 
ever, there is, to our knowledge, no argument for a similar 
effect in the structure functions. 

Assuming the exact result Ca = 1 allows to obtain 
rather accurate values of C.p by using the property of ESS 
(extended self-similarity) |0 plotting Sp as function of 
S3 . (It should be noted that ESS usually results in almost 
perfect scaling behavior, but the scaling coefficients thus 
derived vary in time, hence time averaging is required.) 
The results are shown in Fig. ^, which gives the ESS 
results of Cj" for run 4 (diamonds) and the ESS values 
^+ = Cj^/Ca"^ from ||l| for 2D MHD turbulence (trian- 
gles). (In 2D MHD the absolute values of (p are found 
in [y to decrease with Rm, while the relative values ^3p 
appear to be independent of Rm.) The results indicate 
that in 3D MHD turbulence is less intermittent than 2D, 
but it is more intermittent than hydrodynamic turbu- 
lence, the continuous curve, which gives the She-Leveque 
result Cp^. As shown by Politano and Pouquet |^, the 
She-Leveque concept contains effectively three parame- 
ters: g related to the scaling Szi ~ l^^^, x related to 
the energy transfer time at the smallest, the dissipative 
scales ti ^ l^ , and C, the codimension of the dissipative 
structures. 
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Our results for the 3D MHD case suggest Kolmogorov 
scaling 5 = 3, x = 2/3, while different from hydrody- 
namic turbulence the dissipative structures are sheet-like, 
hence the codimension is C = 1, 



(^f'^^^p/^ + l-il/if/^ 



(8) 



This is the dashed curve in Fig. H, which fits the numer- 
ical values very well. 

In conclusion we have studied the spatial scaling prop- 
erties of 3D MHD turbulence using direct numerical sim- 
ulations with resolution of 512'^ modes. The results indi- 
cate that the turbulence does not follow the Iroshnikov- 
Kraichnan (IK) phenomenology. The energy spectrum 
is consistent with a k^^'^ law. For hyperresistivity the 



structure functions exhibit a clear scaling range yielding 
absolute values of the scaling exponents ^p, in particu- 
lar C3 ~ 1, consistent with recent analytical predictions. 
The scaling exponents agree well with a modified She- 
Leveque model C^^^ = p/9 -I- 1 — (1/3)^/^, correspond- 
ing to Kolmogorov scaling, but sheet-like geometry of 
the dissipative structures. The results are also consistent 
with observations of turbulence in the solar wind, which 
typically show a k~^'' spectrum. 
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TABLE I. Summary of the simulation runs. 



run No 


iV 


v 


Tiv 


Rmx 




1 


128 


1 


3 X lO^-' 


19 


A 


2 


256 


1 


10"^ 


40 


+ 

N 


3 


512 


1 


3 X 10-* 


81 




4 


512 


2 


3 X 10-*^ 


- 





1 


.(a) 


/1-4 
/ 1.2 


Slope 


- 


0.10 


\ 




: / 


1.0 


_ --y~->=^^^^^-^ : 


; 




/ 


0.8 


\ 


- 


0.0 ! 


/ 


0.6 
0.4 


\ - 


: 




100 



10 



00 



000 






1.00 



0.10 ^ 



0.01 




0.01 



0.10 

k 



1.00 



FIG. 1. Scatter plot of the normalized angle-integrated en- 
ergy spectrum compensated with k^'^ from run 3. The dashed 
line indicates the k~^'^ spectrum. 
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FIG. 3. Log-log plot of the time-averaged normalized struc- 
ture functions from run 4. (a) S^{1), (b) S^{1). The inserts 
give the derivatives, the horizontal dashed lines provide most 
probable values of the scaling exponents. 
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FIG. 2. Time-averaged normalized energy spectra compen- 
sated with k^'^ from runs 1-3. 



FIG. 4. Scaling exponents ^jj" for 3D MHD turbulence (dia- 
monds) and relative exponents Cp /Ct for 2D MHD turbulence 
(triangles). The continuous curve is the She-Leveque model 
Cp'", the dashed curve the modified model Cp'^^. 



